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ABSTRACT 


The  Hall  coefficient  and  conductivity  have  been  measured  dur¬ 
ing  pile  irradiation  fora  number  of  zone -refined  p-type  silicon 
crystals  with  initial  resistivities  of  1,  8,  and  100  ohm-cm.  To 
supply  the  magnetic  field  (500  oersteds)for  the  Hall  measurements 
in  the  reactor,  a  small  electromagnet  was  used.  The  conductivity 
of  zone-refined  silicon  shows  much  faster  changes  with  irradia¬ 
tion  than  pulled  silicon  samples  of  equivalent  resistivity.  The 
100-ohm-cm  samples  exhibit  a  monotonic  nonlinear  decrease  of 
in  ar  (conductivity)  vs  <pt  (integrated  fast  flux),  whereas  the  other 
samples  with  initial  Fermi  levels  closer  to  the  valence  band  have 
one  or  two  regions  of  linear  decrease  in  in  a-  vs  <pf  before  the 
nonlinear  decrease  region  is  observed.  The  Hall  mobility  for  the 
100-ohm-cm  samples  decreases  and  becomes  negative  as  a  result 
of  the  carrier  density  decreasing  with  irradiation.  In  the  case  of 
the  8-ohm-cm  sample,  the  Hall  mobility  decreases  with  irradia¬ 
tion,  whereas  the  1 -ohm-cm  sample  shows  no  change  in  Hall  mo¬ 
bility  with  irradiation  up  to  the  maximum  integrated  flux  used  in 
the  present  experiment.  The  origins  of  the  dependence  of  In  a 
on  <P,  as  well  as  the  behavior  of  the  Hall  coefficient  and  Hall  mo¬ 
bility  with  irradiation  are  discussed. 


PROBLEM  STATUS 

This  is  an  interim  report  on  the  problem; 
work  is  continuing. 
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IN-PILE  HALL  COEFFICIENT  AND  CONDUCTIVITY  MEASUREMENTS 
ON  ZONE-REFINED  P-TYPE  SILICON 


INTRODUCTION 

Radiation  damage  studies  of  single-crystal  silicon  have  been  in  progress  for  a  num¬ 
ber  of  years  and  have  revealed  many  interesting  features.  These  include  the  existence 
and  location  of  discrete  energy  levels  (1-10)  formed  in  the  forbidden  band  by  various 
kinds  of  incident  irradiation  (electrons,  deuterons,  neutrons,  and  gamma  rays).  An 
analysis  of  the  radiation  damage  is  usually  accomplished  by  observing  the  postirradia¬ 
tion  temperature  dependence  of  the  Hall  coefficient,  the  electrical  conductivity,  and  the 
optical  absorption.  But  perhaps  the  most  powerful  tool  used  to  study  defect  centers  is 
electron  spin  resonance  taken  in  conjunction  with  optical  absorption.  Using  this  tech¬ 
nique,  it  has  been  shown  quite  conclusively  (11,12)  that  the  Si-A  center,  which  is  respon¬ 
sible  for  the  electron  trapping  level  (0.17  ev  below  the  conduction  band  in  electron- 
irradiated  pulled  silicon),  is  a  combination  of  a  vacancy  and  an  oxygen  atom. 

In  contrast  to  most  of  the  previous  work  which  has  been  carried  out  on  pulled  sili¬ 
con,  having  about  1018  oxygen  atoms  per  cm3,  the  present  experiments  were  set  up  to 
measure  the  effect  of  neutron  irradiation  on  the  Hall  coefficient  and  conductivity  of  zone- 
refined  silicon,  which  has  less  than  1016  oxygen  atoms  per  cm3  (13).  Furthermore,  the 
experiments  were  designed  to  be  carried  out  during  pile  irradiation  in  the  reactor, 
whereas  all  previous  Hall  coefficient  measurements  on  neutron-irradiated  silicon 
(5,6,9,14,15)  have  been  made  after  irradiation.*  Buras  and  Suwalski  (16),  however,  have 
made  in-pile  Hall  coefficient  measurements  on  germanium. 

In  this  report,  after  the  experimental  techniques  of  measurement  have  been  dis¬ 
cussed,  the  experimental  findings  will  be  divided  into  four  main  headings:  the  approach 
of  the  samples  to  intrinsic  conductivity,  the  variation  of  conductivity  with  integrated  flux, 
the  Hall  coefficient  variation  and  hole  removal  rates,  and,  finally,  the  variation  of  Hall 
mobility  as  a  function  of  integrated  flux. 


EXPERIMENTAL  TECHNIQUES 

Samples  with  dimensions  1  x  0.2  x  0.2  cm  were  cut  from  uncompensated  zone- 
refined  single  crystals  of  p-type  silicon  (doped  with  boron)  purchased  from  Merck  and 
Company.  The  long  axis  of  each  sample  was  in  the  (111)  direction.  The  samples  were 
etched  in  CP-4,  sandblasted  with  an  industrial  airbrasive  unit,  and  then  gold  plated.  A 
jig  which  covered  the  ends  and  four  points  (two  for  conductivity  probes  and  two  for  Hall 
coefficient  probes)  of  the  individual  crystals  was  used  while  excess  gold  plate  was  sand¬ 
blasted  away.  The  crystals  were  placed  in  transite  holders  between  the  poles  of  a  small 
electromagnet.  The  magnetic  field  was  maintained  at  a  rather  low  but  sufficient  value  of 
500  oersteds.  Samples  B-l,  D-8  and  G-lOOl  were  all  irradiated  at  100°C  with  the  follow¬ 
ing  exception:  the  temperature  of  D-8  varied  as  the  ambient  reactor  temperature  and 
was  79°C  initially,  93°C  at  90x  1012  n/cm2,  and  100°C  after  195  xlO12  n/cm2.  Samples 


*This  is  true  with  the  possible  exception  of  Wertheim's  study  (5)  of  fission  plate- 
irradiated  silicon. 

'The  number  in  the  sample  designation  refers  to  the  preirradiation  room  temperature 
resistivity  in  units  of  ohiri-cm. 
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H-100  and  J-100  were  irradiated  at  72°C  in  the  presence  of  a  dry  helium  atmosphere. 
Other  details  of  the  irradiation  procedure  are  described  in  an  earlier  paper  (17). 

The  measuring  equipment  consisted  of  a  sensitive  potentiometer  with  an  electronic 
null  detector.  By  means  of  a  suitable  switching  arrangement  the  null  detector  was  con¬ 
nected  either  across  the  galvanometer  output  of  the  potentiometer  or  in  series  with  the 
source.  Very  high  sensitivity  was  obtained  with  the  null  detector  in  series  with  the 
source  and  this  was  especially  needed  as  irradiation  proceeded  and  the  source  resist¬ 
ance  increased.  With  a  100-ohm-cm  sample,  for  instance,  the  source  resistance  varied 
during  irradiation  from  500  ohms  to  above  50,000  ohms.  With  this  apparatus,  five  quan¬ 
tities  were  measured:  the  Hall  potential,  the  potential  across  the  resistivity  probes,  the 
sample  current,  the  temperature,  and  the  current  through  the  magnet  coils.  From  these 
observations,  the  Hall  coefficient  and  resistivity  were  calculated. 

The  flux  for  samples  B-l  and  D-8,  irradiated  at  the  centerline  of  the  core,  was  3.3 
X  10  10  n/cm2  sec  and  was  measured  using  activity  determinations  of  monitors  of  sulfur 
(2.9-Mev  threshold)  pellets  irradiated  with  the  samples.  Samples  G-,  H-,  and  J-100 
were  placed  1  foot,  3.4  feet,  and  4  feet  above  the  core  centerline.  For  these  samples, 
the  flux  was  determined  by  a  separate  experiment,  taking  the  centerline  core  flux  to  be 
that  as  found  for  samples  B-l  and  D-8.  Several  sulfur  pellets  were  placed  at  intervals 
above  the  core  centerline,  and  the  resultant  log  activity  vs  distance  plot  gave  fluxes  for 
G-,  H-,  and  J-100  of  1.1  x  10 10,  1.1  x  109,  and  3.6  x  108  n/cm2  sec,  respectively.*  The 
reactor-power  level  for  samples  B-l,  D-8,  and  G-100  was  100  kw  whereas  for  H-  and 
J-100,  the  power  was  1  Mw.  This  way  of  determining  the  flux  values  for  H-  and  J-100 
assumes  that  the  log  activity  vs  distance  plots  for  the  reactor  operating  at  100  kw  and  1 
Mw  have  the  same  slopes. 


RESULTS  AND  DISCUSSION 
Approach  to  Intrinsic  Conductivity 

The  conductivity  of  all  the  samples  decreases  as  a  function  of  integrated  fast  flux 
and  tends  toward  some  limiting  value  as  shown  in  Fig.  1.  This  value,  o-L  ,  is  tabulated  in 
Table  1  along  with  the  calculated  value  of  cr.i ,  the  intrinsic  conductivity  at  the  tempera¬ 
ture  of  measurement.  Values  of  cr.  were  obtained  from  (18)' 


ai  =  Pi^O,,  +  Mh)  (1) 

where  p;  is  the  density  of  holes  in  the  intrinsic  region  and  //e  and  p.h  are  the  electron 
and  hole  conductivity  mobilities  respectively.  If  one  uses  the  empirical  expressions  for 
the  conductivity  mobilities  and  for  p;  as  obtained  by  Morin  and  Maita  (19),  one  obtains 


Cr 


1.5  x  10 3 3  T3  e 


-  1 . 2  l/kT\  i/2 


(2) 


The  temperature  range  over  which  Eq.  (2)  is  valid  is  that  for  which  lattice  scattering  is 
the  main  scattering  mechanism.  According  to  Conwell  (20)  the  impurity  contribution  to 
the  mobilities  of  the  samples  used  (10 ls ,  1015,  1014  impurity  atoms/cm3)  is  small.  The 
influence  of  radiation-induced  defects  on  the  scattering  may  be  estimated  as  follows.  If 
one  takes  ten  displaced  atoms/cm3  per  neutron/cm 2, *  the  total  density  of  displaced 


*The  assumption  is  made  that  the  integrated  flux,  and  not  the  flux,  is,  the  important  parameter, 
t A  value  of  37  cm- 1  was  calculated  for  silicon  using  the  method  of  Ref.  21.  This  value 
was  arbitrarily  reduced  to  10  cm"1  on  the  basis  of  a  statement  also  in  Ref.  21,  p.  36,  to 
the  effect  that  the  ratio  of  calculated  to  observed  numbers  of  defects  "seems  to  lie  in 
the  range  of  3  to  10  for  most  materials." 
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Fig.  I  -  Conductivity  o-  as  a  function  of  inte¬ 
grated  fast  flux  4>f  for  different  p-type  silicon 
samples  irradiated  at  100°C  (B-l,  D-8,  G-100) 
and  at  72”C  (J-100,  H-100).  The  crosses  are 
values  of  calculated  from  Eq.  (2).  To  ob¬ 
tain  the  abscissas  of  points  in  n./cm2,  multiply 
the  numbers  along  the  abscissa  axis  by  2  x  1014 
(B-l,  D-8),  5  x  1013  (G-100),  2  x  1012  (J-100), 
and  10  12  (H-100).  Some  of  the  points  are  very 
low  integrated  flux  for  sample  B-l  hhve  been 
omitted  for  the  sake  of  clarity. 


atoms  Nd  at  the  end  of  irradiation  is  found  and  the  values  are  given  in  Table  1.  The 
largest  value  of  Nd  is  1.8  x  1016/cm3,  which  is  about  the  same  as  the  initial  chemical 
impurity  concentration  for  B-l.  Since  an  initial  impurity  concentration  of  10  Is/cm3  does 
not  have  much  effect  on  the  mobility,  neither  should  the  radiation-induced  defects,  and 
the  use  of  Eq.  (2)  seems  quite  reasonable.  This  is  more  a  plausibility  argument  than  a 
proof,  since  one  should  really  measure  m-c  and  p.h  as  a  function  of  irradiation. 

It  is  clear  from  comparison  of  crL  and  o-.  in  Table  1  that  samples  D-8  and  G-100 
have  become  intrinsic,  H-100  and  J-100  are  almost  intrinsic,  and  sample  B-l  is  not  in¬ 
trinsic.  It  is  to  be  noted  that  when  the  sample  becomes  intrinsic  the  density  of  displaced 
atoms  is  10  times  the  initial  hole  density  p0 ,  -also  given  in  Table  1.  The  significance,  if 
any,  of  this  factor  of  10  is  not  known  at  present.  The  purer  the  sample,  the  smaller  is 
the  initial  carrier  density,  and  the  smaller  is  the  integrated  fast  flux  necessary  to  make 
the  sample  go  intrinsic.  (Compare  D-8  with  G-100  both  irradiated  at  100°C.)  Also,  the 
lower  the  temperature,  the  smaller  is  the  integrated  flux  needed  to  make  the  sample  be¬ 
come  intrinsic.  (Compare  G-100  with  H-  or  J-100.) 
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To  compare  the  data  of  this  paper  with  that  of  other  workers,  it  is  convenient  to  use 
the  value  of  the  integrated  flux  <pf,  at  which  o-  equals  2°-..;  these  values  are  tabulated  in 
Table  1.  Crawford  and  Cleland  (22)  observed  the  behavior  of  the  approach  to  intrinsic 
conductivity  of  pile -irradiated  n-type  and  p-type  pulled  silicon  having  initial  resistivities 
of  0.8  ohm-cm  and  5  ohni-cm  respectively.  Their  p-type  silicon  irradiated  at  60°C  was 
within  a  factor  of  2  of  being  intrinsic  after  lx  10  17  n/cm2.  This  is  a  factor  of  10 3 
greater  integrated  flux  than  the  value  for  the  8-ohm-cm  sample  irradiated  at  a  slightly 
higher  temperature  in  the  present  experiment.  Longo  (14)  irradiated  2-ohm-cm  n-type 
and  12-ohm-cm  p-type  pulled  silicon  in  a  reactor  and  found  the  samples  to  have  resistiv¬ 
ities  corresponding  to  intrinsic  material  following  the  irradiation  at  45°C.  His  quoted 
integrated  flux  after  irradiation  was  7  X1017  n/cm2.  Thus,  with  respect  to  the  approach 
to  intrinsic  conductivity,  zone-refined  and  pulled  silicon  crystals  of  equivalent  resis¬ 
tivity  behave  in  the  same  manner  except  that  zone-refined  silicon  approaches  intrinsic 
conductivity  with  much  less  irradiation  than  pulled  silicon.  This  agrees  with  a  footnote 
statement  of  Watkins  et  al.  (23),  who  made  measurements  of  electron-irradiated  silicon. 


Detailed  Variation  of  Conductivity  with  Integrated  Fast  Flux 

Sample  B-l  -  Figure  1  shows  that  the  initial  behavior  of  B-l  is  different  from  that 
of  D-8  and  the  100-ohm-cm  samples.  For  sample  B-l,  there  are  initially  two  linear 
regions  of  in  a  vs  <pr  (To  make  certain  that  a  similar  behavior  was  not  missed  in  the 
two  100-ohm-cm  samples  because  of  the  rapid  change  of  a  with  4>f,  a  third  100-ohm-cm 
sample,  J-100,  was  irradiated  at  a  very  low  flux  of  3.6  X  108  n/cm2  sec.  The  conductiv¬ 
ity  data  shown  for  this  sample  was  taken  over  a  period  of  780  minutes,  and  a  monotonic 
nonlinear  decrease  of  In  <r  with  4>,  was  observed.) 

To  help  explain  the  behavior  of  the  conductivity  as  a  function  of  integrated  fast  flux 
for  sample  B-l  (and  for  the  other  samples  as  well),  we  will  examine  the  variation  of  the 
Fermi  energy  with  integrated  flux.  Analysis  of  the  Hall  data  reveals  for  sample  B-l  that 
In  p  vs  <j>f ,  where  p  is  the  hole  density,  also  has  two  linear  regions.  Furthermore,  the 
hole  density  is  related  to  the  Fermi  energy  Ef  by  use  of  (24) 


kT  In 


(,  3/  2 

1 

h2  J  P 


(3) 


where  mh  is  the  density -of-states  mass  for  the  holes  in  the  valence  band  and  Ev  is  the 
energy  of  the  top  of  the  valence  band.  Therefore,  since  In  p  vs  <£f  is  linear,  Ef  vs  4>t 
will  also  be  linear.  In  addition,  we  assume  that  the  lower  half  of  the  forbidden  band  has 
defect  levels  similar  to  those  in  the  neutron-irradiated  pulled  silicon  observed  by  Klein 
and  Straub  (6);  that  is,  there  are  levels  at  0.16  ev  and  0.30  ev  above  the  valence  band.* 
The  density  of  holes  trapped  by  the  defects  at  a  given  temperature  depends  on  the  proba¬ 
bility  of  occupancy  of  the  level  and  the  net  rate  at  which  the  defects  are  formed.  The 
probability  of  occupancy,  shown  in  Table  1,  is  calculated  from 

[. ,  ,<-r->/‘T]'‘ 


and  is  0.17  and  0.94  for  the  0.16-ev  and  0.30-ev  levels  respectively  for  sample  B-l. 
Thus,  if  both  levels  were  formed  at  the  same  rate,  the  0.30-ev  level  would  be  populated 
much  faster  with  holes  than  the  0.16-ev  level,  because  its  energy  is  lower  (as  far  as 
holes  are  concerned)  and  consequently  its  probability  of  occupancy  is  higher.  But  the 


*The  following  discussion  in  the  test  does  not  depend  on  there  being  a  hole  trap  at  0.06 
ev  (6)  above  the  valence  band  since,  the  probability  of  occupancy  of  this  level  is  0.01  for 
B-l,  0.001  for  D-8,  and  much  lower  for  the  100-ohm-cm  samples. 
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0.30-ev  level  is  formed  at  a  lower  rate  (as  will  be  discussed  in  the  following  section) 
than  the  0.16-ev  level.  In  the  case  of  the  1-ohm-cm  sample,  then,  the  Fermi  level  in¬ 
creases  linearly  with  q>,  until  the  combination  of  the  two  factors  above  favors  the  filling 
of  the  0.30-ev  levels  instead  of  the  0.16-ev  levels.  At  this  point,  Ef  begins  to  change 
with  <t> f  at  a  different  rate  and  the  second  linear  region  of  In  o-  vs  <p(  begins.  Through¬ 
out  the  entire  irradiation,  of  course,  Ef  is  increasing  toward  the  center  of  the  gap. 
Finally,  as  saturation  occurs,  there  is  established  an  equilibrium  which  involves  the 
density  of  holes  left  in  the  valence  band  and  the  density  of  the  0.30-ev  levels  being 
formed.  Longo  (14)  observed  more  than  one  slope  in  In  a  vs  integrated  flux  for  0.05- 
ohm-cm-initial-resistivity  (presumably  pulled)  p-type  silicon  irradiated  with  deuterons. 
The  end  of  the  first  linear  portion  of  In  a  vs  integrated  dueteron  flux  in  Longo's  experi¬ 
ment  occurred  when  6  x  10 17  holes/cm3  had  been  removed  from  the  valence  band.  This 
would  occur  when  approximately  10 18  levels/cm3  which  involve  oxygen  atoms*  were 
filled.  In  the  present  work  with  zone-refined  silicon  (oxygen  concentration  less  than 
1016/cm3),  the  end  of  the  first  linear  portion  occurs  when  approximately  lx  10 16  holes/ 
cm3  are  removed  from  the  valence  band.  Thus,  the  density  of  holes  removed  by  the  time 
the  second  linear  portion  of  In  a  vs  <pf  begins  in  both  pulled  and  zone-refined  silicon 
seems  to  be  very  nearly  the  same  as  the  available  oxygen  density.1' 

Sample  D-8  -  For  sample  D-8,  the  Fermi  level  is  initially  0.27  ev  (as  calculated 
from  Eq.  (3)  and  shown  in  Table  1)  above  the  valence  band,  and  one  sees  from  Table  1 
that  the  probability  of  occupancy  of  the  0.16-ev  and  0.30-ev  levels  is  0.026  and  0.72 
respectively.  The  probability  of  occupancy  of  the  0.30-ev  level  is  therefore  28  times 
higher  than  that  of  the  0.16-ev  level,  and  the  difference  in  production  rate  of  the  0.16-ev 
and  0.30-ev  levels  is  not  sufficient  to  counteract  the  much  higher  probability  of  occu¬ 
pancy  of  the  0.30-ev  levels.  Thus,  the  single  slope  in  the  plot  of  In  p  vs  <Pf  or  of  Ef  vs 
<Pf  for  D-8  is  due  to  a  filling  of  the  0,30-ev  levels.  As  saturation  takes  place,  the  change 
of  in  p  with  <£,  is  governed  by  the  rate  at  which  equilibrium  is  established  between  the 
density  of  defects  being  produced  at  the  0.30-ev  level  and  the  density  of  holes  in  the 
valence  band. 

Samples  G-,  H-,  and  J-100  -  In  the  case  of  samples  G-,  H-,  and  J-100,  the  Fermi 
level  is  initially  0.34  ev  or  more  above  the  valence  band  and  the  probability  of  occupancy 
of  the  0.16-ev  defect  level  is  0.0014  for  G-100  and  0.0023.  for  H-  and  J-100.  One  can 
note  that  the  probability  of  occupancy  of  the  0.30-ev  level  is  initially  0.10  or  higher  and, 
although  low,  is  close  to  a  factor  of  100  times  that  of  the  0.16-ev  level.  Thus  for  these 
samples,  as  for  sample  D-8,  the  larger  occupation  probability  of  the  0.30-ev  level  favors 
its  being  filled  in  preference  to  the  0.16-ev  level.  The  fact  that  one  does  not  initially  get 
a  linear  decrease  in  the  curve  of  Ef  vs  <p{  or  in  in  p  vs  <pf  for  these  samples  may  be 


*The  lower  lying  0.16-ev  level'  could  result  from  a  hole  trap  due  to  various  impurities 
(including  boron)  other  than  oxygen.  However,  since  the  conductivity  of  zone-refined 
silicon  changes  much  more  rapidly  with  irradiation  than  that  of  pulled  silicon  of  equiva¬ 
lent  resistivity,  it  seems  reasonable  that  oxygen  is  the  responsible  agent,  simply  be¬ 
cause  this  is  the  primary  difference  between  zone-refined  and  pulled  silicon.  In  addi¬ 
tion,  oxygen  was  chosen  as  the  atom  responsible  for  the  trap  at  0.16  ev  because  of  the 
similarity  between  the  defect  levels  in  the  upper  and  lower  halves  of  the  energy  gap  (8) 
and  because  the  level  0.16  ev  below  the  conduction  band  in  electron-irradiated  pulled 
n-type  silicon  is  due  to  an  oxygen  atom  and  vacancy  combination  (11).  There  is  some 
controversy  as  to  whether  neutron-irradiated  n-type  silicon  has  any  discrete  trapping 
levels.  According  to  Rupprecht  and  Klein  (9)  there  is  a  discrete  level  at  0.15  ev  be¬ 
low  the  conduction  band  but  'Wertheim  (5)  and  Sonder  (15)  claim  there  is  a  continuous 
set  of  levels  in  the  upper  half  of  the  band. 

^The  density  of  holes  removed  during  the  first  linear  portion  is  also  the  same  as  the  ini¬ 
tial  boron  impurity  concentration,  and  the  assumed  hole  trapping  level  at  0.16  ev  above 
the  valence  band  could  be  due  as  much  to  a  boron-involved  level  as  to  an  oxygen- 
involved  level.  See  reference  29  as  to  why  oxygen  was  chosen  as  the  impurity  causing 
the  0.16  ev  level. 
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due  to  a  competition  for  the  holes  between  traps  formed  at  0.30  ev  and  traps  at  higher 
energies. 

If  this  interpretation  is  correct,  one  should  find  no  occupied  0.16-ev  levels  in  the 
1-ohm-cm  sample  after  the  end  of  the  first  linear  portion  of  Ef  vs  4>(1  and  no  occupied 
0.16-ev  levels  at  any  irradiation  dosage  in  the  8-  and  100-ohm-cm  samples.  Two 
recent  papers  by  Vavilov  and  Plotnikov  (25)  and  Plotnikov  et  al.  (26)  partially  confirm 
this  "prediction."  They  made  photoconductivity  measurements  at  100°K  and  concluded 
that  there  are  levels  in  100-ohm-cm  zone-refined  pile -irradiated  (10 13  n/em2  with  an 
average  energy  of  1  Mev)  p-type  silicon  at  0.30,  0.38,  and  0.45  ev  above  the  valence 
band,  and  no  levels  at  0.16  ev.  To  further  confirm  these  ideas,  a  simple  experiment 
could  be  done  to  locate  the  hole  trapping  levels  in  1-ohm-cm  zone-refined  silicon  after 
a  short  enough  neutron  irradiation  such  that  In  <y  is  still  in  the  first  linear  region. 

The  fact  that  the  hole  density  and  conductivity  decreases  immediately  with  irradia¬ 
tion  in  the  case  of  the  100-ohm-cm  samples  may  indicate  that  the  deeper  lying  levels  do 
not  result  from  clusters  of  vacancies  or  interstitials.  Or,  if  the  levels  are  caused  by 
clusters,  the  clusters  are  immediately  formed  when  irradiation  is  begun. 


Hall  Coefficient  Variation  with  Integrated  Fast  Flux; 

Hole  Removal  Rates 

The  variation  of  the  Hall  coefficient  with  integrated  fast  flux,  as  shown  in  Fig.  2, 
can,  of  course,  be  explained  on  the  basis  of  carrier  removal  and  by  resorting  to  the  two- 
carrier  model' as  the  material  becomes  intrinsic. 

The  most  interesting  feature  of  the  Hall  coefficient  data  is  that  R,  initially  positive, 
reaches  a  maximum  as  a  function  of  integrated  flux,  and  then  goes  negative  for  samples 
G-,  H-,  and  J-100.  Similar  results  were  observed  with  pulled  silicon  by  Longo  (14) 
after  much  higher  irradiation  (1017  neutrons/cm2).  With  dueteron  irradiation  of  0.16- 
ohm-cm  p-type  (presumably  pulled)  silicon,  Longo  also  observed  R  to  go  negative  after 
7  x  10M  deuterons/cm2  (9.6  Mev)  at  room  temperature.  Klein  and  Straub  (6),  on  the 
other  hand,  irradiated  60-ohm -cm  (presumably  pulled)  p-type  silicon  in  the  Brookhaven 
reactor  at  dry  ice.  temperature  and  observed  the  Hall  coefficient  and  conductivity  after 
room  temperature  storage  extending  from  several  days  to  a  year.  They  saw  no  change 
in  the  sign  of  R  with  an  integrated  flux  of  10 14  n/cm2.  In  the  present  work  R  becomes 
negative  after  1.7  x  10 13  and  1.4  x  10 13  n/cm2  for  samples  H-  and  J-100  respectively 
and  after  7.2  x  1013  n/cm2  for  sample  G-100.*  These  considerations  seem  to  suggest 
that  R  goes  through  a  maximum  and  becomes  negative  more  rapidly  for  floating  zone- 
refined  silicon  than  for  pulled  silicon  of  equivalent  resistivity. 

Samples  H-  and  J-100,  irradiated  at  different  flux  rates,  show  similar  dependence  of 
R  on  integrated  flux  as  seen  in  Fig.  2.  Sample  J-100  has  a  smaller  maximum  than  H-100, 
probably  because  of  a  small  initial  increase  of  10° C  in  temperature  of  the  former.  After 
a  few  minutes  of  irradiation,  the  temperature  of  J-100  was  back  to  72°C.  However,  the 
positions  of  the  maxima  of  R  occurred  at  the  same  value  of  <£f  and  the  positions  of  R  = 

0  occurred  at  approximately  the  same  value  of  <£f  for  both  samples.  Sample  G-100  was 
irradiated  at  a  different  flux  as  well  as  at  a  higher  temperature  (100°C)  than  either  H- 
or  J-100.  The  maximum  in  R  for  G-100  occurred  at  an  integrated  flux  of  4  times  that 
for  samples  H-  and  J-100.  In  addition  to  this,  R  for  G-100  did  not  become  as  negative 
as  that  of  H-  and  J-100  with  the  integrated  flux  obtained  during  the  experiment.  The 
small  difference  (on  an  absolute  scale)  in  irradiation  temperature  of  30°C  between  G-100 


*Note  from  the  figure  legend  that  relative  to  the  abscissa  scale  for  samples  H-  and  J-100 
the  scale  for  sample  G-100  is  compressed  by  a  factor  of  5. 
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Fig.  2  -  Hall  coefficient  R  as  a  function  of  inte¬ 
grated  fast  flux  <pj  for  different  p-type  silicon 
samples  irradiated  at  100°C  (B-l,  D-8,  G-100), 
and  at  72 ‘C  { J- 100,  H-  100).  To  obtain  the  abscis¬ 
sas  of  points  in  n/cm1,  multiply  the  numbers 
along  the  horizontal  axis  by  2  x  1014  (B-l,.  D-8), 
5  x  lO13  (G-100),  and  1013  (H-100,  J- 100).  Some 
of  the  points  at  very  low  integrated  flux  have  been, 
omitted  for  the  sake  of  clarity. 


and  H-  or  J-100  apparently  made  a  large  change  in  annealing  rate  and  hence  in  the  net 
rate  of  formation  of  the  defect  levels.  This  change  was  sufficient  to  increase  the  re¬ 
quired  flux  necessary  to  bring  the  Hall  coefficient  to  zero  for  sample  G-100. 

Sample  D-8  exhibits  a  much  broader  peak  in  R  vs  4>f  and  remains  positive  during 
the  entire  irradiation  time  used  for  this  sample.  In  addition,  the  other  low  resistivity 
sample,  B-l,  never  exhibited  a  maximum  in  R.  Longo  (14)  observed  r  of  0.16-ohm-cm 
material  to  go  negative  after  a  long  deuteron  irradiation.  This  would  suggest  that  R  for 
the  1-  and  8-ohm-cm  samples  eventually  would  have  gone  negative  if  longer  irradiations 
had  been  used. 

To  examine  the  carrier  removal  rates  we  consider  the  Hall  coefficient  from  which 
the  carrier  density  for  p-type  silicon  is  calculated  using  R  =  r/pe,  where  r  is  the  ratio 
of  yuH,  the  Hall  mobility,  to  ytc,  the  conductivity  mobility.  This  expression  can  be  used 
only  during  the  initial  period  of  irradiation  and  is  not  valid  when  the  sample  is  approach¬ 
ing  intrinsic  behavior.  According  to  the  graph  of  Morin  and  Maita  (19)  for  r  vs  T,  r  is 
about  0.71  and  0.68  for  T  =  72°C  and  100°C  respectively.  Because  of  the  spread  in  their 
data,  a  value  of  0.70,  independent  of  irradiation,  has  been  taken  to  compute  p,  the  hole 
density. 
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The  carrier  removal  rates  are  listed  in  Table  1.  The  values  are  obtained  from  the 
first  few  points  of  a  linear  plot  of  p  vs  <£f,  which  is  consistent  with  the  errors  involved. 
These  values  are,  of  course,  very  approximate  because  of  the  uncertainty  in  the  inte¬ 
grated  fast  neutron  flux.  The  relative  values,  however,  are  significant.  The  initial  re¬ 
moval  rate  for  B-l,  80  cm'  *,  is  due  primarily  to  holes  being  trapped  by  the  0.16-ev 
level.  At  the  end  of  the  first  linear  portion  of  In  o-  vs  the  removal  rate  for  sample 
B-l  has  dropped  to  40  cm' 1  because  the  holes  are  now  preferentially  trapped  by  defects 
with  energy  levels  at  0.30  ev.  In  the  case  of  sample  D-8  the  initial  removal  rate  of  60 
cm' 1  is  lower  than  the  initial  removal  rate  of  B-l.  The  fact  that  the  holes  initially  are 
being  trapped  primarily  by  defects  with  energy  levels  at  0.16  ev  for  sample  B-l  and  at 
0.30  ev  for  sample  D-8  suggests  that  the  0.30-ev  level  defects  are  being  formed  at  a 
lower  rate  than  the  0,16-ev  levels.  The  difference  between  the  initial  removal  rate  of 
D-8  and  the  removal  rate  of  B-l  at  the  beginning  of  the  second  linear  region,  when,  pre¬ 
sumably,  for  both  samples  the  0.30-ev  levels  are  being  populated,  may  be  due  to  the 
lower  initial  temperature  of  D-8.  Further  evidence  for  the  lower  rate  of  formation  of 
the  0.30-ev  levels  lies  in  the  small  initial  removal  rate,  20  cm-1,  for  G-100.  In  this 
sample,  only  the  defects  with  energies  at  0.30  ev  or  higher  can  trap  holes,  as  discussed 
in  the  previous  section.  For  sample  G-100,  the  initial  removal  rate  is  lower  than  for 
D-8  and  for  the  removal  rate  at  the  beginning  of  the  second  linear  portion  of  B-l.  This 
may  be  due  to  the  lower  probability  of  occupancy  (0.1)  of  the  0.30-ev  level  for  G-100  as 
well  as  to  a  competition  for  holes  from  the  valence  band  between  the  0.30-ev-level 
defects  and  higher  level  defects.  This  seems  to  indicate  that  the  higher  level  defects  are 
formed  at  an  even  lower  rate  than  the  0.30-ev-level  defects. 

The  large  difference  in  removal.rates  between  sample  G-100  (20  cm'1)  irradiated 
at  100°C  and  samples  H-100  (115  cm'1)  and  J-100  (140  cm'1)  irradiated  at  72-C  indi¬ 
cates  either  a  large  amount  of  annealing  or  perhaps  a  reduction  in  the  clustering  of  de¬ 
fects  which  may  lower  their  effectiveness  in  trapping.  Additional  work  has  to  be  done  to 
clarify  this  point. 


Hall  Mobility 

The  Hall  mobility  nH  is  shown  in  Fig.  3  as  a  function  of  integrated  fast  flux.  The 
dashed  line  in  Fig.  3  connects  the  100°C  preirradiation  value  of  for  sample  D-8  with 
the  value  observed  after  the  sample  reached  100°C  during  the  irradiation.  The  Hall  mo¬ 
bility  curve  of  H-100  was  very  similar  to  that  of  J-100  and  has  been  omitted  for  the  sake 
of  clarity.  One  sees,  at  first,  that  the  maximum  in  R  as  a  function  of  4>{  does  not  result 
in  a  similar  maximum  in  nH-  The  reason  for  this  is  that  the  increase  in  R  is  compen¬ 
sated  by  the  decrease  in  &  in  the  early  period  of  irradiation.  appears  to  remain  con¬ 
stant  with  4>f  for  sample  B-l,  except  for  the  scatter  in  the  data.  The  Hall  mobility  for 
samples  D-8  and  G-,  H-,  and  J-100  decreases  with  <#>f  and,  for  the  100-ohm-cm  sam¬ 
ples,  goes  through  zero  and  then  becomes  negative.  One  might  at  first  think  that  the 
radiation-induced  defects  may  increase  the  scattering  of  the  holes  and  hence  decrease 
p.n.  However,  as  indicated  above,  if  the  radiation  defects  exert  the  same  influence  on  the 
conductivity  as  substitutional  chemical  impurities,  the  mobilities  should  not  change  much 
with  the  integrated  flux  received  by  the  samples  in  the  present  work.  Also,  since  nH  for 
sample  B-l  remained  constant,  another  cause  for  the  large  decrease  in  nH  must  be 
found. 

As  irradiation  proceeds,  the  hole  traps  cause  a  decrease  in  p.  Consider  the  two- 
carrier  expression  for  the  Hall  mobility: 


=  Rcr 


2,  2 
P  "  nfMe//ih 

r  - 

P  +  n  ne//2h 


(4) 
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Fig.  3  -  Hall  mobility  /xH  has  a  function  of  inte¬ 
grated  fast  flux  4>,  for  different  p-type  silicon 
samples  irradiated  at  100°C  ( B - 1 ,  D-8,  G-100) 
and  at  72  ’C  ( J-  1 00).  To  obtain  the  abscissas  of 
points  in  n/cm2,  multiply  the  numbers  along 
the  horizontal  axis  by  2  x  10  14  (B-l,  D-8),  5  x 
10  13  (G-100),  and  1013  ( J- 100).  Some  of  the 
points  at  very  low  integrated  flux  have  been 
omitted  for  the  sake  of  clarity. 


where  n  is  the  electron  density.  As -p  decreases  and  approaches  n^c2//rh2,  /rH  de¬ 
creases,  goes  through  zero,  and  becomes  negative.  If  the  1-ohm-cm  sample  had  been 
irradiated  long  enough  to  approach  intrinsic  conductivity  with  its  Hall  coefficient  becom¬ 
ing  negative  (both  conditions  were  true  for  the  8-  and  100-ohm-cm  samples),  then  for 
it  would  have  decreased  with  irradiation. 


SUMMARY 

The  variations  in  the  Hall  coefficient  and  the  conductivity  of  zone-refined  silicon 
undergoing  neutron  irradiation  are  similar  to  those  of  pulled  silicon  except  that  the 
changes  with  integrated  flux  occur  much  more  rapdily  for  the  zone-refined  material. 

The  logarithm  of  the  conductivity  shows  a  monotonic  nonlinear  decrease  with  irradiation 
for  the  100-ohm-cm  samples,  a  single  region  of  linear  decrease  for  the  8-ohm-cm  sam¬ 
ple,  and  two  regions  of  linear  decrease  for  the  1-ohm-cm  sample.  For  the  latter  sam¬ 
ple,  this  effect  is  attributed  partly  to  the  filling  of  the  levels  which  may  involve  oxygen 
atoms  in  some  way.  The  initial  hole  removal  rates  are  80,  60,  and  20  cm"1  for  the  1-, 
8-,  and  100-ohm-cm  samples  respectively  at  an  irradiation  temperature  of  100° C,  and 
115  cm" 1  for  one  of  the  100-ohm-cm  samples  at  an  irradiation  temperature  of  72°C. 

The  Hall  coefficient  goes  through  a  maximum  and  decreases  for  the  8-  and  100-ohm-cm 
samples,  becoming  negative  for  the  100-ohm-cm  samples.  In  the  case  of  the  1-ohm-cm 
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sample,  the  Hall  coefficient  increases  and  then  levels  off.  The  Hall  mobility  decrease  in 
the  8-  and  100-ohm-cm  samples  is  ascribed  to  a  decrease  in  the  carrier  density.  The 
Hall  mobility  is  constant  with  irradiation  for  the  1-ohm-cm  sample. 
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